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August 4, 2010 

Abstract. We report a detailed investigation on the generation of pulse pairs during the readout of a 
coherence grating stored in a cold atomic ensemble. The pulse shapes and the split of the retrieved energy 
between the two pulses are studied as a function of the relative intensities of the two reading fields, and 
a minimum is observed for the total retrieved energy. We introduce a simplified analytical theory for 
the process, considering a three-level atomic system, which explains all the most striking experimental 
features. 

PACS. 42.50.Gy Effects of atomic coherence on propagation, absorption, and amplification of light; elec- 
tromagnetically induced transparency and absorption - 32.80.Qk Coherent control of atomic interactions 
with photons 



1 Introduction 

In recent years the techniques for generation, storage, and 
manipulation of light pulses by atomic systems have devel- 
oped rapidly in connection to the interest in constructing 
more general quantum networks for various applications 
in quantum and classical information and metrology [JJ. 
The idea is to combine the photon's capability for fast 
exchange of information with the long lived atomic mem- 
ories in order to allow information to be processed in a 
distributed manner among various sites. In this field, the 
studies of light storage based on the phenomenon of Elec- 
tromagnetic Induced Transparency (EIT) plays a central 
role [2j|3l|4]. Soon after the theoretical proposal by Fleis- 
chhauer et al [5] for such light storage in an atomic ensem- 
ble, the demonstration of this phenomenon was reported 
by Liu et al [5] and Phillips et al [7J. Since then these 
demonstrations have triggered a number of related works 
reporting different aspects of the light storage process [8j 

Parallel to the development of the storage of light pulses 
by atomic ensembles came various proposals for using such 
effect in quantum information applications. The control- 
lable storage of single photons, for example, could be a 
great step in the implementation of linear optics quantum 
computation [13] . The use of atomic ensembles to gen- 
erate and store single photons could also result in large 
efficiency gains for the distribution of entanglement over 
long distances, as in the DLCZ protocol for quantum re- 
peaters |14] . Presently, prove-of-principle experiments have 
already demonstrated both single photon storage 15,16 
and various building blocks of the DLCZ protocol [T7J 
ITS] . Even though such experiments still present efficien- 
cies and coherence times considerably lower than the ones 
required for the actual implementation of the above men- 



tioned quantum algorithms, they demonstrate the versa- 
tility and power, for various applications, of the process 
of light storage by atomic ensembles. 

Our present work explores new possibilities of such 
light storage process involving the extraction of the stored 
information in different directions and the possibility to 
distribute the information in more than one optical mode 
at once. In our scheme, an ensemble of cold atoms is 
initially prepared in a coherent superposition of Zeeman 
ground states by a pair of writing beams with a small angle 
between their respective propagation directions. After the 
turn off of these beams a coherence grating is left in the 
ensemble, storing information related to various aspects of 
the optical state of the writing fields [T21IT9] . In the present 
work, we only use laser pulses as writing fields. However, 
in order to understand the implications of our work, it 
is important to have in mind that this configuration for 
writing optical information in an atomic ensemble is the 
same as the ones employed for both single-photon [TB] and 
squeezed- vacuum [30] storage in cold atoms. 

We emphasize here that the readout of the stored in- 
formation can be manipulated in a quite general way, by 
both simple control of the direction of the retrieved pulse 
or the distribution of the stored information among differ- 
ent optical modes. In order to investigate such processes 
we extract the stored coherence grating using an inde- 
pendent pair of counter-propagating reading beams. In 
this way, by controlling the relative intensity of the two 
read fields we can choose between the extraction of the 
coherence grating as a single optical pulse in one par- 
ticular direction, as a single optical pulse in the oppo- 
site direction, or as a pair of counter-propagating light 
pulses of arbitrary relative amplitudes. Such pulse pair 
generation from a stored coherence grating was first re- 
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ported in Ref . . Here we approach this process under a 
more general framework, investigating the transition be- 
tween the more well known processes of Delayed-Four- 
Wave-Mixing [T21I2T] and EIT [51[7] readouts and such 
pulse pair generation. We also introduce a simple, ana- 
lytical theory that accounts for all major features of this 
double-field readout, and use it to discuss various aspects 
of the process. 

In the following, Section [5] introduces the general idea 
of our scheme and the theory to model it, together with 
the analysis for various aspects of the pulse-pair genera- 
tion process. Section [3] describes our experiment and its 
main results. It also provides the comparison between ex- 
perimental results and the theory of Sec. [2] Such compar- 
ison reveals a quite satisfactory agreement, particularly if 
we keep in mind the relative simplicity of the theoretical 
model. Finally, in Sec. 2] we provide a summary of our 
results and the perspective for future works. 



2 Theory 

As mentioned above, we consider here an ensemble of cold 
atoms sequentially excited by two pairs of laser fields, see 
Fig. 1. The first pair is composed by the two writing fields 
W and W', and the second pair by the two reading fields 
R and R'. The atoms have a lambda configuration for 
the relevant level structure, with two degenerate ground 
states, | la) and |16), and one excited state |2). These states 
correspond to specific Zeeman sublevels of two hyperfine 
states, and thus the optical fields connecting them have 
specific circular polarizations. The polarization of fields 
W and R', connecting states \la) and |2), is then a + . 
Fields W and R, on the other hand, have a~ polarization 
and connect levels |16) and |2). These polarizations are 
defined in the reference frame of the atom. Fields W and 
W' have an angle 9 between their respective propagation 
directions. The angle is small enough for us to neglect any 
relative polarization change of the beams due to this angle. 
Fields R and R' have opposite directions, but propagate 
in the same line as W. 

The pair of writing fields act on the atomic ensemble 
for a long time before they are turned off, leaving a coher- 
ence grating printed in the ground state of the ensemble 
(Fig. la). This first stage of coherent preparation of the 
system is discussed in Sec. 12.11 After some storage time 
i s , the two read fields are turned on, resulting in the emis- 
sion of two other fields, D and D' (Fig. lb). These two new 
fields are generated as pulse pairs with relative intensities 
given by the relative powers of the two reading fields. They 
carry the information about the coherence grating origi- 
nally stored in the atomic ensemble. This reading process 
will be discussed in detail in Sec. 12.21 



2.1 Coherent Preparation 



and prepare the necessary notation for the following sec- 
tions. In this way, we consider fields W and W' as plane 
waves propagating in the directions specified by the wave- 
vectors kw = kwz and kw', respectively, with electric 
fields given by 



E w = S w e l(kwz ' uwt) a + 
E w > = £we i ^'- r - w ^' t h 



(la) 
(lb) 



where £w and £w> represent the amplitude of each field, 
and uw and u>w' their frequencies. As discussed above, 
vectors kw and k\y' form a small angle 6. 



(a) Writing 




(b) Reading 



D' 




Oscilloscope 




Fig. 1. (colors online) Schematic representation of the (a) 
writing and (b) reading processes. Top: propagation directions 
of the writing (W,W), reading (R,R'), and diffracted (D,D') 
fields. Base: level structure considered in the theory, with indi- 
cation of the energy levels connected by each field. 7^2 and 7 
are the relaxation rates of the excited state and of the ground- 
state coherence, respectively. 



We are interested in the situation where the rate 7 of 
decoherence between the ground states is much smaller 
then the spontaneous decay rate Z22 of the excited state, 
and both fields W and W are exactly on resonance. In 
this case, the steady-state coherence p\ a lb between levels 
I la) and \lb) is given by 



Pla.lb 



ir 



\n 



(2) 



with 



id 2i i a £we tkv 



id 2 ,ib£we 



(3a) 
(3b) 



The coherent preparation process that we consider here 
was introduced and discussed in detail in Ref. [T2] . In this 
section, then, we will simply review these previous results 



the Rabi frequencies related to fields W and W, respec- 
tively, and dij the dipole moment between levels i and 
j. In the same situation, the populations of the ground 
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states are given by 

Pla,la 
PlbAb 



fi- 



ll " 



\(2 W \ 2 + \O w > 
\n w ,\ 2 



\Q 



+ \Q 



w 



(4a) 
(4b) 



If we wait a time t s > > 1 /I22 after the fields W and 
W' are turned off, the only non-zero elements of the den- 
sity matrix p s specifying the state of the system at t s are 



P\a,la — 2 2^P^aAa P\h,\b) 1 
P\b,\b = Q _ ~2^Pla,la ~ Plb,lb) > 



PlaAb — Pla,lb e 



(5a) 

(5b) 
(5c) 



This gives then the state of the atomic ensemble, with the 
stored coherence grating, right before the readout fields 
are turned on at an arbitrary time t s » 1//22- 



2.2 Pulse pair generation 

The reading process starts once fields R and R' are turned 
on. As discussed above, these two fields propagate in op- 
posite directions along the z axis and their electric fields 
are given by 



E R = £ R e*- k *'- u * t '>d-- , 
E R , =£ R ,e i ^' z - u ^a + , 



(6a) 
(6b) 



with £ R and Sri their constant amplitudes, ui R and ui R i 
the laser frequencies, and k R and k R / specifying their 
wave- vectors along the z axis. We have then k R / = kw = 
—k R . Defining the Rabi frequencies fin and fl R i analo- 
gously as in Eqs. (J3jt) and {3p), we can then write the 
following Bloch equations for the time evolution of the 
system in the readout process 



dpi a ,la 
dt 



dpu, 



lb 



dt 



dcria,: 
dt 



da 1 



b:2 



dt 



= [-Q R > (Tla,2 + C.C.] 

H — PlaSa — Plb,lb) j 

= [-QR(Tlb,2. + C.C.] 

^ — Pla.la ^ Plb.lb) , 

= —f} R ,(l — 2p la ,j a — Plb,lb) + & R pla,lb 

2~ CT la,2 ! 

= — — ^Plb.lb — Pla.la) + ^ R 'Plb,la 

^22 
7r a lb,2 , 



dpia. 



lb 



dt 



— CT2,16 — ^R Olaa 



(7a) 

(7b) 

(7c) 

(7d) 
(7e) 



where we neglected the decoherence rate 7 of the ground 
state, considering a fast readout with !7 R , Q r i » 7. 
The laser fields are again at exact resonance, and o\ a ,2 = 
Pia^e - " 75 '* and uib.2 = Pib,2 e ~" Rt are the slowly varying 
coherences. 

In order to obtain the temporal shape of the retrieved 
pulses, we need to analytically solve this system of eight 
equations for both <J\ a ^{t) and 01^2 (i), considering the 
initial conditions as the stored state discussed previously, 
i.e., Eqs. {Sh.) to {5b) plus oi a , 2 (0) = oi hj2 (0) = 0. An 
important step for such analytical solution is the definition 
of the following new variables: 



and 



Q = ^R&lb,e , 
P = ^R'Vla,2 , 
T = f}R>f2 R Pla,lb 



Qr = Q + Q* , 

Qi ' Q Q 1 

p r =p+p* , 
Pi=p-p* , 

T r = T + T* , 
Ti = T-T* . 



(8a) 
(8b) 
(8c) 



(9a) 
(9b) 
(9c) 
(9d) 
(9e) 
(9f) 



In terms of this last set of variables, Eqs. {7^.) to ([7b- ) split 
in two decoupled systems, which we call just "System 1" : 



195 

Ti - -^P, 



dt 

dQ 1 = _ T £22 
dt 2 ^ ' 

^§ = \n R ,\ 2 Q t -\n R \ 2 p t ; 

dt 

and "System 2" : 

dpiaM r> 1 P22 1 -. v 
-T, = ~Pr + —(1 - PlaM - PXb,Xb) , 



(10a) 
(10b) 
(10c) 

(11a) 



dp 



16,16 



"~2~~(1 — Pla,la " Plb,lb) , 



dt 

= -2|J? fl/ | 2 (l - 2p laM ~ Plb,lb) + T r —P r > 



(lib) 

[ Pr , 

(lie) 

— - = — 2\Sl R \ 2 (l — 2plb,lb — Pla,la) + T r 2~ < ^ r ' 

(lid) 

(He) 



dT r 
~dT 



= -\Q R ,\ 2 Q r -\Q R \ 2 P r 



The solution of System 1 is quite straightforward. Con- 
sidering the initial conditions and defining 



I t = I R + I R > = 



&{\Q R \ 2 + \f2 R ,\ 2 ) 



r 2 
1 22 



(12) 
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an adimensional quantity roughly proportional to the total 
intensity of light exciting the atoms, we obtain 



with 



fr(t) 



p i {t) = T i {Q)f T {i)/r 22 , 

Qi(t) = -Ti(0)f r (t)/r 22 , 



-r 22t /4 sinh (yi _ 2I t r 22 t/4) 



VT=2T t /4 



(13a) 
(13b) 



(14) 



Note that I R and I R i give the respective intensity of the 
fields in units of the saturation intensity of the transitions 
as defined in Ref. [22] . In this way, It is the total intensity 
of the field in units of the respective saturation intensities. 

The solution of System 2 is more elaborated [53] , and 
we have to define yet another adimensional quantity 



Id — Ir — Ir/ 



(15) 



proportional now to the difference between the intensities 
of fields R and R' . Taking into account the initial condi- 
tions, we have then 



P r (t) = r 22 f r (t) 
+ r 22 g r (t) 



. 4i* 
SI t 



'{PlaAa ~ PlbAb) + 



It - Id(Aa,la - PlbAb) + 16 



r P (o) 



r 2 

1 22 

(16a) 



Ut) = -r 22 f r (t) 



\Pla,la Plb,lb) ' 



it r 22 



r229r(t) 



Sit 



It - Id{p e laM - Pi M6 ) + 16^^ 



22 

(16b) 



with 



9r(t) 



2r 3 + 2ri + 1 
2(r 2 - ri)(r 2 - r 3 ) 
2r 2 + 2ri + 1 



(e 



riT 22 t _ „r 2 T 22 t\ 



e r\r 22 t _ e r 3 r 22 t\ 



2(r 3 - ri)(rs - r 2 ) 
The coeficients r%, r 2 , and r% are given by 



(17) 



1 

2 ' 



r 2 = — - 



T\ = s + v 
(s + v) 1 



-~2 + — (S ~ W)Z ' 



(s + v) 1 V3 



(s - v)i , 



with 




In the same way as f r (t), g r (t) depends then only on I t 
and -T 2 2- 

With Pi, P r , Qi, and Q r , it is then possible to invert 
Eqs. © and © and obtain directly 



Clo,2(t) 
0-16,2 (*) 



P r (t)+Pj(t) 

Q r ft)+Q,ft) 



(18a) 
(18b) 



which gives then the medium polarization responsible for 
the emission of the D, D' pulse pair. 



2.2.1 Equal populations in the ground states 

An important situation to be analyzed is when \ fiw\ = 
which leads to p\ a \ a = p\h\b- This case corre- 
sponds to a maximum visibility of the stored coherence 
grating, and to simpler expressions for both eri a;2 and 

<T la , 2 (t) = -i9r(t)^e- ik *'- r 
^1 22 

~i l 9 ' [fr{t)I R +9r{t)Iw] 
Altl 22 

x e -i(k n + k w —k w /)-r 

, j R ,\n R \e-^° rj M 

+ * ^r-^ [Jr{t) - g r {t)\ 



<rib,2(t) 



2i t r 22 

x — i(2h R /— k R — kw+k w /)-r 
<^ 22 



(19a) 



[fr{t)Iw + 9r(t)h 



2i t r 22 

x e -i{k R ,-k w +k w ,yr 

■ 1 TprT, Urit) - g r (t)\ 

£Hi 22 

x — i(2k R — k R i+kw— k w ,)-r 



(19b) 



where we used fix = i\flx\e lkx " r , with X being R 7 R', W, 

A 



or W. We also substituted pi a sb(0) — Pi a u, = e _7 *< 
as is the case once | flw \ = \ Qw 1 \ ■ 

As discussed in more detail in Ref. [12], the second 
term in Eq. flT9"k ) is the one responsible for the D field 
back propagating with respect to W . Correspondingly, 
the D' field comes from the second term in Eq. (|19b ). 
which generates a field propagating in the direction of W . 
The first term on these two equations correspond to stim- 
ulated emission on the respective transitions, contributing 
to fields on the directions of R and R' . The third term, 
for the geometrical configuration of the beams we con- 
sider here, is not phase-matched and does not result in a 
propagating field. 
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For the sake of the following discussions, we only need 
to consider then the electric field amplitudes of D and D': 

I ftp I 

En(t) « ±rj?- [f r (t)I R + g r (t)lR'} , (20a) 

22 

Ed- (t) ex [f r (t)I R > + 9r(t)I R ] . (20b) 

22 

For detectors with a fast response, the respective mea- 
sured signals [12 should be given by 



S D (t) =A D I -§ \f r (t)I R + g r (t)I R ,\ 2 , 



(21a) 



Iw 



S D 'it) = A D ,-2 r \f r (t)I R , +g r (t)I R \ z , (21b) 

with Ajj and Ar>' two proportionality constants. Equa- 
tions (l2~Th ) and (|21b ) will be directly compared to the 
experimental results of Sec. [3] and theirs applicability will 
be discussed there in more detail. Examples of theoretical 
pulse shapes can be found in Fig. A physical situation 
where only the f r (t) function determines the temporal dy- 
namics was explored in Ref. |12j . 



2.2.2 Energy in the diffracted modes 

Another important quantity to be investigated is the re- 
trieved energy in modes D and D', since these are directly 
related to the retrieval efficiency in the reading process for 
the information stored in the atomic ensemble [23] ■ They 
will be given by [12] 



U D 


= / S D {t)dt 
Jo 

= B D T 4 f°|/r(t)^ + 

h Jo 








g r (t)i R '\ 2 dt , 


(22a) 




= / S D >(t)dt 
Jo 

t rco 
= B D'^T \fr(t)I R . 

1 t Jo 








+ g r (t)I R \ 2 dt , 


(22b) 



with Bo and Bp* two new proportionality constants. In 
this way, the total energy retrieved will be 



U T (t) = U D + Ujy 



(23) 



All these quantities (Ud, Ud>, and Ut) will be experimen- 
tally investigated in the next section. For the following 
calculations and analysis, we are going to consider com- 
pletely symmetrical conditions for the generation of D and 
D', so that B D = B D ,. 

In general, Ut depends only on the three intensities 
It, I R , and I R ' . The way we chose to explore this quantity 
was to fix the total intensity and change only the relative 
ratio between the two reading beams, i.e., we studied Ut 
as a function of I R in a situation where I R / — I t — I R . 
We obtain then graphs as shown in Fig. [2Ji, where we 
plotted Ut as a function of I R for five different values of 



It, going from 0.01 to 100. These plots clearly demonstrate 
that the total diffracted energy presents a minimum at 
I R = I R > = once I t gets large enough (I t > 1). This 
behavior indicates then a specific kind of saturation in this 
system, which prevents the whole retrieved energy to be 
evenly distributed among these two modes for large total 
intensities. 




0.0 0.2 



Fig. 2. (a) Normalized total energy in the diffracted modes 
D,D' as a function of the R beam intensity Ir, in a situation 
where the R' beam intensity is I R i = It — Ir. It is then kept 
constant for each curve, with its value taken, from top to bot- 
ton, as 0.01, 0.1, 1.0, 10, and 100, respectively. Note that the 
minimum always occurs at Ir — I R i — It/ 2. The curves are 
normalized by their value at Ir — 0. (b) Total energy in the 
diffracted modes at the minimum as a function of the total in- 
tensity It- For each It, the curve is normalized by the Ut value 
with Ir = and Iri — I t . 



From the theoretical expressions deduced above, this 
behavior comes as a consequence of the fact that f r (t) ~ 
g r (t) in the limit I t << 1. If we make f r — g r in Eqs. (j22"|) . 
it is straightforward to see that Ut becomes independent 
of the relative ratio between I R and I R i and, consequently, 
does not present a minimum anymore. On the other hand, 
in the limit It > 1 the function |g r (i)| 2 becomes consis- 
tently smaller than \f r (t)\ 2 , resulting in the overall depen- 
dence of Ut with I R , I R i . The minimum formation as It is 
increased can be better visualized in Fig. [2Jo, in which the 
solid line represents the minimum Ut(Ir = I R ' = It/2) 
normalized by the value of the total diffracted energy when 
all the reading power is in just one reading beam, i.e., 
normalized by Ut{I R — 0, I R i — It) as in Fig. [2ji. We see 
then that we go from no minimum at It << 1 to a value 
asymptotically approaching 0.277 at I t >> 1. The inter- 
pretation of this specific value for the asymptotic limit of 
the minimum is still under investigation. 

At first sight, the dependence of the total diffracted en- 
ergy with the relative intensities of the reading beams may 
be attributed to any process that can affect the diffrac- 
tion efficiency in a similar way, and which is maximized 
when the two reading intensities are equal. One such pro- 
cess could be spontaneous emission, which might erase the 
stored information, therefore reducing the conversion effi- 
ciency of energy of the beams R,R into the beams D,D' . 
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However, although the excited state population depends 
locally on the relative reading beam intensities, one can 
easily show (Appendix [A} that the total amount of spon- 
taneous emission in the atomic ensemble depends only on 
the total intensity I t and cannot be responsible for the 
observed dependence. Similarly, a calculation of the to- 
tal energy emitted by the atomic ensemble in the form 
of stimulated emission also leads to an expression which 
depends only on I t (see also Appendix [A} , not presenting 
a minimum at Ir = Ir> = -ft/2. Therefore, we are left 
with the coherent process associated with the third terms 
of Eqs. (fT5)) . These terms represent a spatially depen- 
dent polarization induced simultaneously by both read- 
ing beams and are responsible by a portion of the lo- 
cal atomic excitation. Although they represent a process 
which is not phase matched for the reading beam config- 
uration of Fig. [1] and, therefore, do not lead to the gen- 
eration of a coherent beam, their excitation leads to a 
decrease in the conversion efficiency of the reading beams 
into the diffracted D,D' modes. Furthermore, it can also 
be shown (Appendix [B]) that the total irradiated energy 
associated with these terms scales with the product of 
the two reading beam intensities, reaching a maximum 
when Ir = Ir> = It/2. Thus, we attribute the predicted 
minimum as being due to excitation of these non-phase- 
matched optical coherences in the medium. 

3 Experiment 



quarter- waveplate so to acquire opposite circular polariza- 
tions and are incident in the MOT forming a small angle 
(9=lCT 3 rad. 



(a) 




AOM-W 




Detector 2 



t \J2 



Diode Laser 



AOM-R 



Pulse Generator 



(b) 



W, W" "ON" 


"OFF" 








R, R' "OFF" 1^4 — 


"ON" 




( 


) time 



For the experiment we employ a sample of cold cesium 
atoms obtained from a magneto-optical trap (MOT), with 
the trapping beams detuned about 12 MHz from the cy- 
cling transition 6S 1/2 (F = 4) -> 6P 3/2 (F' = 5), and the 
recycling repumping beam resonant with the open transi- 
tion 6S 1/2 {F = 3) -> 6P 3/2 (F' = 3). Initially the atoms 
are prepared in the lower hyperfine ground state by switch- 
ing off the repumping beam for about 1 ms and waiting for 
an optical pumping period (induced by the non-resonant 
trapping beams) of about 50 fj,s after which most of the 
atoms are pumped into the &S\/2(F = 3) ground state. 
During this optical pumping period the MOT quadrupole 
magnetic field is also switched off. We use three pairs of or- 
thogonal Helmholtz coils to compensate for residual mag- 
netic fields. The light storage experiment is performed 
using light from an external cavity diode laser (ECDL) 
locked on the cesium closed transition 6Si/2(F = 3) — > 
6P 3/2 (F' = 2). 

The main parts of the experimental setup are describe 
in Fig. [3^,. In order to control the intensities of the incident 
writing beams W and W we use a pair of acoustic-optical 
modulators (AOM-W) operating in opposite diffraction 
orders so to keep the laser frequency resonant with the 
atomic transition. Similarly, another pair of AOMs ( AOM- 
R) allowed for the independent control of the intensities 
of the reading beams R and R' in relation to the writing 
beams. As shown in Fig. [3^,, half-waveplates and polar- 
izing beam splitters allowed us to control of the relative 
intensities of the writing beams as well as of the reading 
beams. The two writing beams W and W pass through a 



Fig. 3. (a) Setup for the simultaneous detection of the D 
and D' pulses. Main components: external cavity diode laser; 
M, mirror; BS, beam splitter; PBS, polarizing beam splitter; 
A/2, half-wave plate; A/4, quarter-wave plate; AOM, acoustic- 
optical modulator. After detection, both signals are sent to a 
four-channel 400 MHz oscilloscope, (b) Time sequence for writ- 
ing (W,W') and reading (R,R') fields. t a is the storage time 
counted from the turn off of W ,W' up to the turn on of R,R' . 
The origin of the timescale is the one used in Figs. [4] and [6] 

The writing pulses are applied for a sufficiently long 
period, approximately 40 fjs, so to allow for the creation 
of the stationary ground-state Zeeman coherence grating. 
The two counter-propagating reading beams R and R 1 are 
switched on after a controllable storage time t s (fixed in 
2 /is in our case) measured from the instant of the turn- 
ing off of the writing beams, Fig. [3}). As can be seen in 
the experimental scheme depicted in Fig. the reading 
beams also have opposite circular polarizations, with the 
R' beam having the same polarization and propagating in 
the same direction as the writing W beam. Under these 
conditions, we were able to retrieve the two counterprop- 
agating Bragg diffracted pulses, labeled D and D', which 
are detected by two photodiodes. It is worth mentioning 
that we have verified that the two retrieved pulses also 
have opposite circular polarizations, consistent with the 
conservation of angular momentum in the whole process. 

In Fig. [4] we show the temporal shape of the retrieved 
pulses for different intensity ratios of the reading beams 
R and R'. These results demonstrate that we can control, 
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in a straightforward manner, the relative intensity of the 
two pulses. However, they cannot be directly compared to 
the theory of the previous section, since we employed two 
relatively slow detectors (PDA36A from Thorlabs) in our 
experiments. Their response time was about 0.5 fis, which 
gives then a lower bound for the measured temporal width 
of the pulses. 
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Fig. 4. Diffracted signals as a function of time for differ- 
ent relations between the intensities of the reading fields: (a) 
Ir'/Ir = 0.02, (b) I R ,/I R = 0.57, (c) I R ,/I R = 1.04, and (d) 
Ir'/Ir = 1-76. The intensity of the writing fields, W and IT', 
is 1.8mW/cm 2 each. The total intensity of the reading fields 
is the same, 4.5mW/cm 2 , for each frame. 



to have the same value for both transitions. In the exper- 
iment, we employed a total intensity of 4.5 mW/cm 2 . The 
theoretical value is consistent with the experimental one, 
as being just above the saturation intensity. Notice, how- 
ever, that a more direct comparison of the two values is 
not possible under our simplified three-level model, since 
we do not take into account the whole Zeeman structure 
of the real atom. 




Fig. 5. Energy of the diffracted pulses, D and D' , as a func- 
tion of the normalized intensity of field R. The squares (circles) 
correspond to the experimental points for D (D'). In the exper- 
iment, the total reading intensity is kept constant and equal 
to 4.5mW/cm 2 . The intensity is 1.8mW/cm 2 for both writ- 
ing fields. The solid curves are the corresponding theoretical 
results from Eq. (I22|l . with It = 1.3. The error bars were ob- 
tained from three independent measurements under the same 
conditions. 



With this in mind, the first experimental results we 
compared to the theory were the ones for the energy of 
the pulses as a function of the relative reading intensity, 
Fig. [5] This energy depends on the integral of the tempo- 
ral shape of the pulses, which we expect to be less sensitive 
to the detector's time response. Figure [5^ plots then the 
energy extracted in each of the diffracted fields, Ud and 
Ud' , as a function of the intensity of field R normalized by 
the total intensity I t . The normalization of the Ut values 
is done by dividing them by the value of Ut with Ir = 
and I hi — I t . For the normalization of the experimental 
values in the plots of Fig. [3J we have carefully accounted 
for the 50% intensity loss in the retrieved D beam caused 
by the extra beam splitter in its pathway as well as for 
all other losses along the paths of both D,D' beams. The 
two writing fields have the same intensity (1.8mW/cm 2 ), 
so we can compare the experimental results to the simpli- 
fied theoretical expressions of Sec. 12.2.11 The theoretical 
curves in Fig. [5] were generated by adjusting I t in order to 
obtain the best fit to the experimental points. We consid- 
ered then I t = 1.3 in units of the saturation intensity of 
the transitions [see Eq. (fl"2"j) and Ref. [22]]: considered here 



The measurements in Fig. [5^ demonstrate that the to- 
tal energy extracted by the reading process presents a min- 
imum when the two reading fields have the same intensity, 
as expected from the theory introduced above. This is pre- 
sented most clearly in Fig. [5b, a plot of the sum of the two 
curves in Fig. [5ji. We observe then a reduction of about 
30% in Ud + Ud 1 when Ir = Igi. The relatively simple 
theory employed here reproduces quite well this behavior, 
indicating that the relevant physical process behind it is 
properly taking into account by the theory. 

As discussed in connection to Fig.[H the measurements 
of the experimental pulse shapes are limited to a minimum 
temporal width of about 0.5 fis due to our detector's re- 
sponse time. In this way, the theoretical pulse shapes for 
fast detectors, as given by Sd and Sd' in Eq. (|2"Tj) . cannot 
be directly compared to the experimental curves. Even 
though, the theory properly describes the split of the ex- 
tracted diffracted fields among the two modes, D and D'. 
In order to illustrate the kind of pulse shapes obtained 
from such theory under the conditions of Fig. [5j we pro- 
vide in Fig.[S]the pulse shapes of D and D' for each of the 
intensity relations of Figs.0^ to[3Ji. 
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Fig. 6. Diffracted signals as a function of time according to 
Eq. (|2ip and for different ratios of the reading intensity: (a) 
Iw/Ir = 0.02, (b) I R ,/I R = 0.57, (c) I R ,/I R = 1.04, and (d) 
Ir'/Ir — 1-76. The total intensity, normalized by Z22, is kept 
constant and equal to It = 1.3. 



4 Conclusions 

In the present work we reported a detailed study of the 
process of pulse pair generation from a coherence grating 
stored in an atomic ensemble. We showed that it is possi- 
ble to generate pulse pairs with arbitrary relative ampli- 
tudes depending on the intensity relation of the two beams 
employed in the reading process, where the stored grating 
is mapped into the two diffracted fields. We also observed 
the formation of a minimum for the total extracted energy 
when the two reading beams have the same intensity. We 
provided a simple, analytical theory that accounted for 
all these features observed in the experiment. From this 
theory, we obtain that such minimum formation is the 
result of saturation of the system by the reading fields, 
which excite other polarization modes that do not lead to 
diffracted, propagating fields. 

From a more general perspective, the process discussed 
above provides a continuous interpolation between two 
more well known reading processes, the ones based on 
Delayed-Four- Wave-Mixing [T2"Il2T] and EIT 0[7], respec- 
tively. This general pulse pair generation process may be 
used as a means, for example, to deterministically gener- 
ate complex, correlated pairs of optical fields. A study of 
the correlation properties of these pulse pairs is currently 
in progress. It may also be employed to distribute the in- 
formation present originally in a single optical field among 
two other fields, after its storage and possible manipula- 
tion in the atomic ensemble. This may be an important 
tool to explore new possibilities in quantum networks. 
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A Spontaneous and Stimulated Emissions 

The total amount of spontaneous emission at a given time 
is given by /22P22) with P22M the total amount of pop- 
ulation in the excited state. In order to obtain an ex- 
pression for it, our first step is the derivation of P22 = 
1 — Pia,ia — Pib,ib for a particular atom. From Eqs. (fTTa ) 
and (fTTb h we have 



Pr ~t~ Qr 



Using Eqs. 
then obtain 

= 



P22 
dt 

under the condition 



A2P22 • 
C2w\ 



(24) 



\n 



w 



Ii 



16T r (0) 



r 2 
1 22 



^22 e 



-r 22 t 



8 



dt' . 

From the definition of T r (t) with the initial state given by 
Eq. (St): 



T r (0) 



cos(fcfl/z + k 



w 



(26) 



When integrating over the whole ensemble [as done to 
go from Eqs. (|T9| to Eqs. f|20[) ] . such term multiplied by 
the cossine of a function of the spatial grating averages 
to zero. The total population in the excited state is then 
proportional to 

PLW oc f* gr (t')e r ^'dt' . (27) 



Notice that this expression depends only on the total in- 
tensity it employed in the reading process and not on the 
relative intensity of R, R' . 

The calculation of the energy emitted by the atomic 
ensemble in the form of stimulated emission is even more 
straightforward, since it corresponds to the first term on 
the right in Eqs. (TH)a ) and (IT5b ). Following similar steps 
as for the derivation of Ud, Ud>, we obtain then expres- 
sions for the energy emitted by the atomic system in the 
reading modes: 



Ur cc I r / \g r (t)Vdt , 
Jo 

poo 

Ub><xI r > / \9r(t)\ 2 dt. 
Jo 



(28a) 
(28b) 



These result in an expression of the total energy emit- 
ted by stimulated emission, 



U$ = U R + U W oc I t \ \g r (t)\ 2 dt 



(29) 



that is also independent of the relative intensity of the 
reading beams. 
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B Third terms of Eqs. GWl 



In order to understand how the third terms of Eqs. ([T 
may affect the energy diffracted in the D,D' modes, we can 
calculate the total energy scattered incoherently by those 
terms of the induced polarization, or even consider a differ- 
ent experimental situation for which these terms are phase 
matched (as in the beams configurations of Refs. [23] and [HDjg 
In both cases, the energy scattered by each of these po- 
larization components are given by: 



17. 



18. 



U a = B, 



U h = B, 



If 



\f r (t)-g r (t)\ 2 dt 



II 



o 

2 poo 



\fr(t) ~ 9r(t)\ 2 dt 



for the transitions starting at level la and 16, respectively. 
Under usual conditions, we should have completely sym- 
metrical emissions for both transitions and, consequently, 
B a — £>;,. The total energy related to these terms would 
be then 



Ut = B, 



IrI 



RJ-R) 



\f r (t)-g r (t)\ 2 dt 



Now we do have a process for which the maximum emitted 
energy would occur for the same conditions at which a 
minimum is observed in the total energy of modes D, D'. 
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